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In this talk, the program package GoSam is presented which can be used for the automated 
calculation of one-loop amplitudes for multi-particle processes. The integrands are generated in 
terms of Feynman diagrams and can be reduced by d-dimensional integrand-level decomposition, 
or tensor reduction, or a combination of both. Through various examples we show that GoS AM 
can produce one-loop amplitudes for both QCD and electroweak theory; model files for theories 
Beyond the Standard Model can be linked as well. 
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1. Introduction 



Recently we have seen tremendous progress in the automation of NLO multi-leg calculations 
[jl|, ^, [|, |J |5| |6] 0| . In addition, public NLO tools containing a collection of hard-coded individual 
processes, like e.g. MCFM [|, g] and VBFNLO [0, [II], ^ have been developed and constantly 
enriched. The matching of partonic NLO predictions to parton shower Monte Carlo programs has 
also seen major advances, see e.g. the developments in Sherpa Jl3|,[l4|], MC@NLO/Herwig++ [15 



lfj, POWHEG-Box []17[ |18p, POWHEL [Jig gOj, |21j]. Therefore we are progressing well towards 
the aim that basically any process which may turn out to be important for the comparison of LHC 
findings to theory can be evaluated at NLO accuracy. 

In this talk, we present the program package GoSam [|I|] which allows the automated gener- 
ation and evaluation of one-loop amplitudes for multi-particle processes. To produce results for a 
certain process specified by the user, there is an "input card" to be edited specifying the details of 
the process. Then the user can launch the generation of the source code and its compilation, with- 
out having to worry about internal details of the code. The individual program tasks are steered via 
python scripts. 

The integrands of the one-loop amplitudes are generated in terms of Feynman diagrams, 
thus allowing to perform symbolic manipulations of the corresponding algebraic expressions prior 
to any numerical step. For the reduction, the program offers the possibility to use either a. d- 
dimensional extension of the OPP method [E2, 23, 24], as implemented in SAMURAI pBj], or tensor 



reduction as implemented in golem95C [|2g, |27|] interfaced through tensorial reconstruction at the 
integrand level [|28j] . 

The program can be used to calculate one-loop corrections within both QCD and electroweak 
theory. Beyond the Standard Model theories can be interfaced using FeynRules []2"9|] or LanHEP 



The Binoth Les Houches interface [ 3 1 ] to programs providing the real radiation contributions is also 
included and has been tested in various examples. 



2. Description of the program 



2.1 Generation of the source code for the amplitude 



For the diagram generation we use the program QGRAF Q32[], supplemented by Python rou- 
tines to further analyse or filter the diagrams, which allows for example to drop diagrams whose 
colour factor turns out to be zero, or to determine the signs for diagrams with Majorana fermions. 
The information about the model is either read from the built-in Standard Model file or is generated 



from a user defined LanHEP D30rj or Universal FeynRules Output (UFO) 029Q file. The program 
also produces a LTjjX file which contains graphical representations of all diagrams together with a 
summary of conventions, e.g. the helicity and colour basis. 

The amplitude is generated in terms of algebraic expressions based on Feynman diagrams and 



then processed with a FORM program, using spinney [33] for the spinor algebra. In GoSam we 
have implemented the 't Hooft-Veltman scheme (HV) and dimensional reduction (DRED). In both 
schemes all external vectors (momenta and polarisation vectors) are kept in four dimensions, while 
internal vectors are kept in the J-dimensional vector space (d = 4 — 2s). For the loop momentum 
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q we introduce the symbol \i 2 = —q 2 , such that 

q 2 =q 2 +q 2 =q 2 -H 2 , (2.1) 

where q lives in 4 dimensions and the (D — 4) -dimensional orthogonal projection is denoted by q. 

To prepare the numerator functions of the one-loop diagrams for their numerical evaluation, 
we separate the symbol /I 2 and dot products involving the momentum q from all other factors. 
All subexpressions which do not depend on either q or \i 2 are substituted by abbreviations, which 
are evaluated only once per phase space point. Each of the two parts is then processed by hag- 



gles [|34[], which generates optimised Fortran95 code for the numerical evaluation. For each 



diagram we generate an interface to Samurai [||], golem95C [|7|] and/or PJFRY [|J5| ||]. 
Our standard choice for the reduction is to use SAMURAI [p5|], which usually provides a fast and 
stable reduction of the amplitude to a set of coefficients of basis integrals in most of the phase 
space. Furthermore, Samurai monitors the quality of the reconstruction of the numerator. In 
GoSam we use this information to trigger an alternative reduction with either golem9 5C [ p7| ] or 



PJFRY [35] whenever these reconstruction tests fail. This combination of on-shell techniques and 
traditional tensor reduction is achieved using tensorial reconstruction at the integrand level 
The tensorial reconstruction not only can cure numerical instabilities, but in some cases also can 
reduce the computational cost of the reduction. Since the reconstructed numerator is typically 
of a form where kinematics and loop momentum dependence are already separated, the use of a 
reconstructed numerator tends to be faster than the original procedure, in particular in cases with a 
large number of legs and low rank. 

2.2 Rational terms 

Terms containing the symbols \i 2 or e in the numerator of the integrands can lead to a so-called 



R2 term [37], which contributes to the rational part of the amplitude. As we start from Feynman 
diagrams, we generate the R2 part along with all other contributions without the need to seperate 
the different parts. In addition, we provide two different ways to calculate R2, an implicit and an 
explicit construction, using the fact that there are two ways of splitting the numerator function: 

^(q,H 2 ,£) = ^6(<?,At 2 ) + e^(-?,At 2 ) + £ 2 ^(<?,M 2 ) (2-2) 
or, alternatively, 

^(q,H 2 ,E) = J^(q)+J^(q,^ 2 ,e). (2.3) 



The implicit construction uses the splitting of Eq. ( |2.2| ) and treats all three numerator functions Jfi 
on equal grounds. Each of the three terms is reduced seperately in a numerical reduction and the 
Laurent series of the three results are added up taking into account the powers of e. 

The explicit construction of R2 is based on the assumption that each term in jV in Eq. ( p3| ) 
contains at least one power of \i 2 or e. The expressions for those integrals are relatively simple 
and known explicitly. Hence, the part of the amplitude which originates from jV is computed 
analytically whereas the purely four-dimensional part jV is passed to the numerical reduction. 

In the user input card, possible options for R2 are r2=implicit, explicit, off and 
only. Using r2=only discards everything but the R2 term and puts GoSam in the position 
of providing R2 terms to supplement other codes which work entirely in four dimensions, provided 
they use the same gauge. 
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2.3 Conventions 

To be specific, we consider the case where the user wants to compute QCD corrections. In the 
case of electroweak corrections, the analogous conventions apply except that the strong coupling 
g s is replaced by e. In the QCD case, the tree-level matrix element squared can be written as 

\^\l ee = ^^o = (g s ) 2h -a , (2.4) 

where b = is also possible. The matrix element at one-loop level, i.e. the interference term 
between tree-level and one-loop, can be written as 



(2.5) 



£ e 2 



A call to the subroutine Samplitude returns an array consisting of the four numbers (ao,co,c_i,c_2) 
in this order. The average over initial state colours and helicities is included in the default setup. 
Renormalisation is included depending on the options chosen by the user, for a more detailed de- 
scription we refer to [[]]]. 

3. Installation and Usage 
3.1 Installation 

The user can download the code GoSAM either as a tar-ball or from the subversion repository 
at http : / /pro ject s . hepf orge . org/ go s am/ . The build process and installation of Go- 
Sam is controlled by Python Distutils, while the build process for the libraries Samurai 
and golem9 5C is controlled by Autotools. To install GoSAM, the user needs to run 

python setup. py install -prefix MYPATH. 

For more details we direct the user to the GoSam reference manual coming with the code. 



On top of a standard Linux environment, the programs FORM [38], version > 3.3, and QGRAF [32] 
need to be installed on the system. Further, at least one of the libraries SAMURAI [25] and go- 
lem9 5C p7| ] needs to be present at compile time of the generated code. For the user's con- 
venience we have prepared a package containing Samurai and golem95C together with the 



integral libraries OneLOop [|39|], QCDLoop Q4Q] and FF Q41J]. The package gosam-contrib- 
1 . . tar . gz containing all these libraries is available for download from 

http :/ /pro ject s . hepf orge . org / go s am/. 

3.2 Usage 

In order to generate the code for a process, the user needs to prepare an input file, called 
process card in the following, which contains 

• process specific information, such as a list of initial and final state particles, their helicities 
(optional) and the order of the coupling constants; 

• scheme specific information and approximations, such as the regularisation and renormali- 
sation schemes, the underlying model, masses and widths which are set to zero, the selection 
of subsets of diagrams, etc; 
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• system specific information, such as paths to programs and libraries or compiler options; 

• optional information for optimisations which control the code generation. 

If the process card is called gosam.in, it can be invoked by the command gosam.py 
go s am. in. All further steps are controlled by the generated make files; in order to generate 
and compile all files relevant for the matrix element one needs to invoke 

make compile 

The generated code can be tested with the program matrix/test. f90. The following sequence 
of commands will compile and run the test program: 

cd matrix 
make test, exe 
. / test . exe 

The numbers printed by the test program are, in this order, a^, cq/oq, c-i/oq, C-2/ao and the pole 



parts calculated from the infrared insertion operator [u2, 43]. One can also generate a pictorial 



representation of all generated diagrams using the command 
make doc 

which generates the file doc/process. ps 

3.3 Interfacing the code 

3.3.1 Using the BLHA Interface 



The so-called Binoth Les Houches Accord (BLHA) [ pl[ ] defines an interface for a standardized 
communication between one-loop programs (OLP) and Monte Carlo (MC) tools. GoSam can act 
as an OLP in the framework of the BLHA, such that the calculation of complete cross sections is 
straightforward. 

In general, the MC writes an order file, called for example olp_order . lh, and invokes the 
script gosam. py as follows: 

gosam.py — olp olp_order . lh 

The invocation of gosam.py generates a set of files which can be compiled with a generated 
make file. The BLHA routines are defined in the Fortran module olp_module but can also be 
accessed from C programs. 

3.3.2 Using External Model Files 



GoSam can also make use of model files generated by either Feynrules [44] in the UFO 



format [ ]29[ ] or by LanHEP [|30|]. The particles can be specified by their PDG code. Details about 
how to import these model files are described in the GoSam reference manual. Precompiled 
MSSM_UFO and MSSM_LHEP files can also be found in the subdirectory examples /model. 
The examples directory also contains examples showing how UFO or LanHEP model files are 
imported. 
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4. Examples and validation 

The code generated by GoSam has been compared to a considerable number of processes 
available in the literature, as listed in Table [j]. Many of these processes are also included as exam- 
ples in the code, including reference values. 




50 100 150 200 250 300 350 400 ' -4 -3 -2 -1 1 2 3 4 

P, (GeV) 



(a) Transverse momentum of the leading jet for W +jet pro- (b) Pseudorapidity of the leading jet for W + jet production at 
duction at LHC with ^ = 7 TeV. LHC with y/s = 1 TeV. 

Figure 1: NLO calculation of W~+jet production at LHC using GoSam interfaced with SHERPA via the 
Binoth Les Houches interface and compared to MCFM. 



As an example for the Binoth Les Houches Accord interface of GoS AM we present results for 
the QCD corrections to W~ + 1 jet production, obtained by linking GoSam with SHERPA [|(]]. 
Results for the transverse momentum and rapidity distribution of the leading jet are shown in 
Figs. 1(a) and 1(b) . The comparison with MCFM ||] shows perfect agreement. Furthermore, 
SHERPA offers the possibility to match NLO calculations with a parton shower [O, [H| ]. 



5. Conclusions 

We have presented the program GoSam which can produce code for one-loop amplitudes for 
multi-particle processes in an automated way. The program is publicly available at 
http : / /pro ject s . hepf orge . org / go s am/ and can be used to calculate one-loop ampli- 
tudes within QCD, electroweak theory, or other models which can be imported via an interface to 
LanHEP or FeynRules. Monte Carlo programs for the real radiation can be easily linked via the 
interface defined by the Binoth Les Houches Accord. 

The amplitudes are generated in terms of Feynman diagrams and can be reduced by unitarity 
based reduction at integrand level or traditional tensor reduction, or a combination of the two 
approaches. Further, the user can choose among different libraries for the master integrals and 
different regularisation and renormalisation schemes. 

The calculation of the rational terms can proceed either together with the same numerical 
reduction as the rest of the amplitude, or before any reduction, using analytic information on the 
integrals which can potentially give rise to a rational part. This feature also allows to use the code 
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p6 , 47], own analytic calculation 



[52] 



[53] 



§1,MCFM[|,|] 



for IR poles: [M ||] 



[57, v3] 



Table 1: Processes for which GoSam has been compared to the literature. 



for the calculation of rational parts only. Moreover, the GoSam generator can produce code for 
processes which include unstable particles, i.e. intermediate states with complex masses. 

GoSam is very well suited for the automated matching of Monte Carlo programs to NLO 
virtual amplitudes, and therefore hopefully will be used widely as a module to produce differential 
cross sections for multi-particle processes which can be compared directly to experiment. 
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